Pulsar Velocities without Neutrino Mass 
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We show that pulsar velocities may arise from anisotropic neutrino emission induced by resonant 
conversions of massless neutrinos in the presence of a strong magnetic field. The main ingredient 
is a small violation of weak universality and neither neutrino masses nor magnetic moments are 
required. 
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One of the most challenging problems in modern as- 
trophysics is to find a consistent explanation for the high 
velocity of pulsars. Observations JlJ show that these ve- 
locities range from zero up to 900 km/s with a mean 
value of 450 ± 50 km/s . An attractive possibility is 
that pulsar motion arises from an asymmetric neutrino 
emission during the supernova (SN) explosion. In fact, 
neutrinos carry more than 99% of the new-born proto- 
neutron star's gravitational binding energy so that even 
a 1% asymmetry in the neutrino emission could generate 
the observed pulsar velocities. To find the origin of such 
asymmetry is, however, not a minor task. One intrigu- 
ing possible explanation to this puzzle may reside in the 
interplay between the parity non-conservation present in 
weak interactions and the strong magnetic fields which 
are expected during a SN explosion. Although several 
possible realizations of this idea in the framework of the 
Standard Model (SM) of particle physics have been al- 
ready explored |^,|| a conclusive solution of the problem 
is still lacking, and there is some motivation for looking 
also at solutions that involve physics beyond the SM. 

Recently, several neutrino conversion mechanisms in 
matter have been invoked as a possible engine for power- 
ing pulsar motion. Kusenko and Segre proposed a mecha- 
nism H based on MS W conversions || . The idea is based 
on the observation that the strong magnetic field present 
during a SN explosion gives rise to some angular depen- 
dence of the matter induced neutrino potentials raj. As 
a consequence, in the presence of non- vanishing v r mass 
and mixing the resonance sphere for the v e — v r conver- 
sions is distorted. If the resonance surface lies between 
the v T and v e neutrino spheres, such a distortion would 
induce a temperature anisotropy in the flux of the escap- 
ing tau-neutrinos produced by the conversions, hence a 
recoil kick of the proto-neutron star. In order to account 
for the observed pulsar velocities the required strength of 
the dipolar component of the magnetic field between the 
two neutrino spheres must exceed 10 15 Gauss M or even 



larger H|. Another crucial ingredient in this mechanism 
is the neutrino squared mass difference, Am 2 > 10 4 eV 2 , 
which leads to m„ T > 100 eV or so, assuming a negligi- 
ble v e mass. This is necessary in order for the resonance 
surface to be located between the two neutrino-spheres. 
It should be noted, however, that such requirement is 
at odds with cosmological bounds on neutrinos masses 
unless the r-neutrino is unstable. 

Akhmedov, Lanza and Sciama [^| proposed a similar 
pulsar acceleration mechanism based on resonant neu- 
trino spin-flavour precession (RSFP) |ic[| . The mag- 
netic field not only affects the medium properties, as 
in the matter density in the MSW case, but also in- 
duces the spin-flavour precession through its coupling 
to the neutrino transition magnetic moment fllTf . The 
lowest magnetic field strength required is B > 2 x 10 16 
Gauss, as long as the neutrino magnetic moment exceeds 

^ > 10~ 15 Ms- 

In this letter we investigate the relevance, for pulsar 
motion, of a different kind of neutrino conversions not 
requiring any neutrino mass nor magnetic moment. The 
basic mechanism was proposed over ten years ago in ref. 
p"2[ . In contrast to the physics of the solar neutrino prob- 
lem, the new mechanism was shown to be potentially 
relevant for supernova physics. This has been recently 
studied in more detail in ref. |p"3|| where stringent limits 
have been derived. 

The simplest underlying particle physics model that 
realizes this new conversion mechanism postulates the 
existence of two new SU(2) <g) (7(1) singlet leptons for 
each generation of leptons, in such a way that lepton 
number symmetry is exact in the Lagrangian |]l4f . These 
extra states can arise in various extensions of the SM, 
such as superstring models However, the model is 

very interesting on its own right, both conceptually as 
well as phenomenologically [^6| . As a result of the pos- 
tulated lepton number symmetry neutrinos remain mass- 
less to all orders of perturbation even after the gauge 
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symmetry breaking. However, unlike the situation in the 
SM there is a non-trivial Kobayashi-Maskawa-like mixing 
in the weak leptonic charged current jL7|. The simplest 
such scheme contains three two-component gauge singlet 
neutral leptons S added to the three right-handed neu- 
trino components v c present in SO(10). For defmiteness 
we consider this model at the SU(2) (8 U(l) level. The 
assumed conservation of lepton number leads to a neu- 
tral mass matrix with the following texture in the basis 
{v, v c 7 S): 



D \ 

D T M 
M T J 



(1) 



where the Dirac matrix D describes the coupling between 
the weak doublet v and the singlet v c , and whereas M 
connects the singlet states v c and S. It is easy to see 
that, as expected, the three conventional neutrinos re- 
main massless, while the other six neutral 2-component 
leptons combine into three heavy Dirac fermions. This 
model offers a viable alternative to the see-saw model. 

The phenomenological implications of this picture are 
manifest when considering the resulting charged- current 
Lagrangian in the massless- neutrino sector: 
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£ C c = -T=Wpe a i J 'Y l iK ai UiL + h.c. , 



(2) 



where a = e, h,t, i = 1,2,3 and the mixing matrix K is 
given as OJlT 



K = RM, 



(3) 



For defmiteness and simplicity we will confine ourselves 
to the case of two neutrinos. In this case R is a 2 x 2 
rotation matrix, 



R 



cos v sin t 
— sin 9 cos ( 



and where the diagonal matrix, 



M- 




(4) 



(5) 



Note that K is not unitary, since it is a sub-matrix 
of the full rectangular matrix including also the heavy 
states The matrix M describes the effective 

non-orthogonality of the two neutrino flavours, i.e., 



{v e \v^ T ) = — sin9 cos9(M 2 ~M£ 3 )- The non-diagonal el- 



ements of the matrix K cannot be rotated away through 
a redefinition of the massless-neutrino fields. In this way 
a non-trivial mixing arises among the massless neutrinos. 

The corresponding form of the neutral-current La- 
grangian for the massless-neutrino sector is fllJTT 



NC 
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where 



P = K^K = M 2 



(6) 



The matrix P is diagonal but generation-dependent, sig- 
nalling the violation of weak universality. It is also con- 
venient to define 



JV? = (1 + A?)" 



1,2(3), 



(8) 



where the hi parameters reflect the deviation from the 
standard neutrino coupling. Since no oscillations be- 
tween two strictly massless neutrinos can develop in vac- 
uum, it follows that laboratory limits on the leptonic 
mixing angle 9 are very weak. However, it will be suffi- 
cient for our purposes to assume that the mixing angle 9 
is very small. In this way we have Vi ~ v a [a = e,/x(r)], 
so that h 2 ~ h 2 . The parameters h 2 are also constrained 
experimentally. There have been extensive experimental 
studies of the constraints on ft, 2 . For the third generation 
one can still allow ft 2 values in the range of a few per- 
cent jl8| , whereas the constraints on h\ and /i 2 are more 
stringent. For this reason we consider from now on the 
case of v e v T conversions. 

In the flavour eigenstate-type basis v a defined in ref. 
0|, the Schroedinger evolution equation describing the 
neutrino propagation in matter can be written as 



dr 



= V2G F ^- 
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where A e ^ T are the neutrino amplitudes, Gf is the Fermi 
constant, p is the matter density, and is the nucleon 
mass. The entries of the evolution Hamiltonian are now 
given by 

h e = y' e {M e c 2 + m t « 2 ) 2 - \y^c 2 + A/ -2 * 2 ), 

H er = [Y^MeC 2 +M T S 2 ) - \Y^{M e +M T )](Mr-Me)sC, (10) 
H T = Yis 2 C 2 {M 2 -M 2 ) - \y: i {M 2 S 2 +M 2 C 2 ) 



Y' = Yn + Y° COS < 



(11) 



where 

Y' e = Y e -y e °cos^ 

where Y e = - n ^ n and the superscript zero refers to 
the contribution from the lowest Landau level |l{|. In 
cq. ( |l0| ) we used the shorthand notation s = sin 9 and 
c = cos 9. The effect of the magnetic field enters through 
the neutrino potentials [H that affect the entries of the 
Hamiltonian. This can be parametrized by defining ef- 
fective electron and neutron fractions Y£ and Y r [ which 
now contain an anisotropic term cos <fi where <fi is the an- 
gle between the neutrino propagation direction and B. 
Using eq. (|l0|) the resonance condition, V e =V T , reads 
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or, in a more transparent form, 



(12) 



F e (l + A e cos^) = 7 ? F n ^l-A e ^cos^ . (13) 

where A e = "n? e /n e where the superscript zero again the 
contribution to the electron fraction coming from the low- 
est Landau level. The parameter 77 defined as 
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(14) 



measures the deviation from lepton universality which, 
for the (is e , v T ) system may reach 77 ~ 10~ 2 or so Jl8| . 

In order to establish that massless neutrino conversions 
can play a role explaining the origin of pulsar velocities, 
we need to verify that the resonance condition eq. ( |l2] ) 
can indeed be fulfilled in a SN environment between the e 
and the r neutrino-spheres. The mean resonance position 
is obtained by averaging eq. ( |l3| ) over cj), giving 



Y e = r,Y n 



(15) 



Apart from small B induced corrections to Y e , the con- 
dition eq. ([l5|) coincides with the free-field resonance 
condition given in rcf. ]l3| which allows us to apply here 
some of the arguments used there. As Y n m 1 in the 
core of the SN, we see that condition in eq. ( p"5| ) can be 
fulfilled for experimentally allowed rj values if Y e < 10~ 2 . 
This is indeed possible close to the neutrino-spheres as a 
consequence of the strong deleptonization taking place in 
that region during the Kelvin-Helmholtz cooling phase. 
A rough estimate of the value of Y e between the neutrino- 
spheres can obtained following ref. |l3| . From the approx- 
imate chemical equilibrium for e~, p, n, and u e , we have 
/i e - +fJ. p ~ where p e - is the electron chemical poten- 
tial and we have set \i Ve ~ 0. Using Boltzmann statistics 
for non-relativistic nucleons and the above chemical equi- 
librium condition , we can write 



Y, 
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exp( Me -/T) + l' 



(16) 



neglecting the neutron-proton mass difference. 

The chemical potential for relativistic and degenerate 
electrons near the neutrinosphere is given by 

Me - « (37T 2 n e ) 1/3 w 51.6(r e />i2) 1/3 MeV, (17) 

where P12 is the matter density in units of 10 12 g cm~ 3 . 
For typical conditions between the two neutrino-spheres, 
T ~ 4 MeV and p 12 ~ 1 + 10 we find y e ~5x 1(T 2 ~ 
5 x 10~ 3 . These analytical results are in good agreement 
with numerical SN models. Therefore, we can expect res- 
onant massless-neutrino conversion to occur between the 
two neutrino-spheres for a range of values of the param- 
eter 77 which is not experimentally excluded. 



The adiabaticity of resonant massless-neutrino conver- 
sions can be easily verified by looking at the probability 
for v e <-» v T and v e <-> v T conversions, given by |L3] 



/ vr 2 Sr 

P = 1 — exp — 

2 L° 



(18) 



1 — exp 



-32pi 2 (ro)( ^=3 



V 1 / h Yc ( r o) 



lcm 



sin 2 26 



Resonant conversion is adiabatic if P ~ 1, hence when- 
ever the resonance width 5r = 2hy c sin 29 is larger than 
the conversion length at resonance L^ m . From eq. (fisl ) we 
see that such a condition is fulfilled, under the conditions 
present between the neutrino-spheres, if sin 2 29 < 1CP 6 . 

The anisotropy in the total momentum of the escaping 
neutrinos can be computed as in [||^] starting from the 
resonance condition eq. (|l3|). Since, as we discussed 
Y e <C Y n in the region of our interest, we neglect the 
second term on the right side of eq. (]l^). Following Q] 
we parametrize the resonance surface equation as follows 



r{4>) = ro + S cos 4> 



(19) 



where ro is the radius of the free-field resonance sphere 
defined by eq. (|l5|). The displacement 5 is found by 
subtracting the resonance condition computed for <f> = tt 
with the same computed in <f) — —ir. We find 

dY 

Y e {r + 5)- Y e (r - S) ~ — -28 = 2\ e Y e {r ) . (20) 

dr 

Defining the Y e variation scale height by 
'dlnFV 1 



hv = 



dr 



we can rewrite eq. (EOp in the more compact form 



5 = Xehyjro) 



(21) 



(22) 



The deformation of the resonance sphere implies an an- 
gle dependence in the temperature of the escaping v T 's, 
hence to an asymmetry in the momentum that these neu- 
trinos carry away from the SN. This is given by B 



A A; 
IT 



1 J* F v cos 



>sm 



sin 
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-Skr, 



(23) 



where the temperature variation scale height hx is de- 
fined in analogy to the definition of hy e given in eq. j2l|). 
The factor 1/6 in eq. ( p^ ) accounts for the fact that, out 
of the six neutrino and anti-neutrino species, we are as- 
suming that only the v r carries a momentum anisotropy. 
We note, however, that an extra factor of two may be 
gained in eq. ( |23] ) if, at the same time, v T emission suf- 
fers resonant conversion between v e and v T spheres. This 
is indeed possible as the resonance condition eq. (|l^) for 
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v e — v T and v e — v T conversions coincide. This is a char- 
acteristic feature of this mechanism Jl2[ which is in sharp 
contrast with the case of MSW conversions || , where ei- 
ther neutrinos or anti- neutrinos can resonantly convert, 
but not both. 

By substituting eq. (p2[) in cq. d2^) we finally get 



Afc 



2 h Y 
9 ~hi 



(24) 



Numerical SN simulations [g0[ typically give ft.y e //it ~ 1 
between the two neutrino-spheres. Hence, we see from 
eq. ( pi[ ) that in order for massless neutrino resonant 
conversions to account for the observed pulsar velocities 
one needs A e ~ 5 x 1CP 2 in the region between the two 
neutrino-spheres. Using the expression for the polariza- 
tion A e given in |pl[| we determine the required value 
of the magnetic fields strength as « 10 15 Gauss, which 
seems reasonable from the astrophysics point of view p2| , 
and lower than required in the MSW [0 and RSFP |||] 
neutrino conversion mechanisms. 

This is illustrated in Fig. 1 for three different choices of 
pY e at resonance consistent with SN models. The three 
lines correspond to pY e values (in g/c.c.) 5 x 10 9 (solid), 
1 x 10 10 (dashed) and 5 x 10 10 (dotted) at resonance. In 
Fig. 1 we assume that both v T and v T carry a momentum 
anisotropy. The kick velocity values given in the right 
ordinate assume a pulsar mass of 1.4 M® and a total 
energy released by all neutrino species of 3 x 10 53 erg. 
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FIG. 1. Magnitude of kick velocity versus magnetic field, 
in units of the critical field B c = 4.4 x 10 13 Gauss. 

In short we have proposed a viable scheme for generat- 
ing pulsar velocities from anisotropic neutrino emission 
induced by resonant conversions of massless neutrinos in 
the presence of a strong magnetic field. Neither neutrino 
masses nor magnetic moments are required and the pa- 
rameters required are totally consistent with cosmology 
(e.g. nucleosynthesis) and astrophysics (e.g. supernova 



physics). However, our proposal rests on the idea that 
there is a small violation of universality in the weak inter- 
action. Again, we are safely within the presently allowed 
region of parameters. However, this mechanism would 
become inconsistent should the sensitivity of weak uni- 
versality tests improve significantly. 
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